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Abstract 

Semi-analytical expressions are suggested for the temperature dependence of those 
combinations of transport coefficients which govern the fission process. This is based 
on experience with numerical calculations within the linear response approach and the 
locally harmonic approximation. A reduced version of the latter is seen to comply with 
Kramers' simplified picture of fission. It is argued that for variable inertia his formula 
has to be generalized, as already required by the need that for overdamped motion the 
inertia must not appear at all. This situation may already occur above T ~ 2 MeV, 
where the rate is determined by the Smoluchowski equation. Consequently, comparison 
with experimental results do not give information on the effective damping rate, as often 
claimed, but on a special combination of local stiffnesses and the friction coefficient 
calculated at the barrier. 

PACS number(s): 24.10.Pa, 24.75.+i,25.70.Jj,25.70.Gh 
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1 Introduction 



It is of considerable interest to understand the temperature dependence of transport prop- 
erties associated with slow collective motion of large scale. Fission is a prime example, and 
indeed, for this case there is growing experimental evidence 0-0 that damping effectively 
increases with T. One often tries to characterize this feature by one parameter, the effec- 
tive damping rate rj which is related to the equation of average motion for a locally defined 
damped oscillator 



through 



V = , ^, , • (2) 

I C I 



The q = Q — Qo measures the deviation of the collective variable Q from some fixed value Qq. 
In the following we will also need other combinations of inertia M, stiffness C and friction 
7, namely 

T-coU = T7T\= — . T^in = — = and W = — — . (3) 

icoU 7 J- kin M 

The Tcou sets the scale for (local) relaxation of collective motion in a given potential of (local) 
stiffness C. The r^in, on the other hand, measures the relaxation of the kinetic energy to 
the equilibrium value of the Maxwell distribution. Typically, for slow collective motion we 
expect this time to be smaller than the former. The limit of overdamped motion applies for 
T'kin ^ TcoU- Using the rj introduced in (^ the following useful relation for their ratio is easily 
verified 



2r] = = ujTcoii = ^ . (4) 

For a positive stiffness (C > 0) and underdamped motion, the zu would be the frequency 
of the vibration and the Fkin its width. It should be noted that in the literature a different 
notation is sometimes used where 7 stands for rj, and often the Fkin//i is referred to as (3. 

To understand the dynamics in phase space one also needs the diffusion coefficients. At 
small temperatures they may deviate from the classic Einstein relation (see |^), but these 
finer details will be neglected here. For such a situation the fission decay rate is commonly 
calculated within Kramers' "high viscosity" limit 0, for which the dependence on friction 
is given by 

"^^'■^^^^ - ^,-vb^(\/^,+vy . (5) 



Rrivb = o) 

Here, the index "b" refers to the fact that the transport coefficients are to be calculated at 
the barrier. It has been reported, see e.g. Fig. 5 of |^, experimental data to suggest, when 
analyzed on the basis of (||), the 77 to be negligibly small at very low temperatures but to 
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rise more or less sharply around T ~ 1 MeV. This result is in qualitative agreement with 
microscopic calculations of the transport coefficients within linear response theory §|, 0, 
although caution is still warranted here. Let us leave aside the fact that at lower temperatures 
the one dimensional potential may attain more structure than that found in one minimum 
and one barrier, which is the picture underlying formula (^) and which for the sake of 
simplicity shall be applied in the sequel. As will be demonstrated below, even then it 
is not permissible to entirely parameterize the truly complicated transport process by the 
single quantity rj. Rather, other combinations of M, 7 and C are needed for more realistic 
descriptions. Moreover, one should be guided by the theoretical fact, that not all transport 
coefficients are equally well accessible both theoretically as well as numerically, as is so for 
the inertia, for instance. 

Unfortunately, when physicists address transport problems, all too often one disregards 
the importance of the inertia, and in particular its variation with the collective coordinate. 
Indeed, in studies based on the Caldeira-Leggett Hamiltonian (see e.g. [l^) or in applications 
of the Random Matrix Theory (see []TT|) the inertia is the one of the unperturbed collective 
part of the total Hamiltonian, treated as an (unknown) parameter. In the case of nuclear 
physics the situation is more complicated. Firstly, there, no unperturbed inertia exists at 
all; it may only show up in the final effective equation of motion as one manifestation of 
the existence of collective dynamics. Secondly, the inertia M may depend sensitively on 
the collective degree of freedom. This feature is already well known from the traditional 
case of undamped motion at zero thermal excitation [jl2|. There is no a priori reason why 
this should be different at finite temperature, with perhaps two exceptions or modifications. 
With increasing T the M gets close to the liquid drop value |jl3|, which only varies smoothly 
with Q and which is quite small. Simultaneously the friction strength increases, so that one 
may quickly reach the situation of overdamped motion, for which no trace of the inertia can 
be seen anymore. Such features have been seen within the linear response approach (see 
0), but to the best of our knowledge no other transport model has so far addressed this 
question. 



2 Rate formulas 

Like in the analysis underlying p, |^ we want to make use of a simple formula for the decay 
rate. In slight modification of Kramers' classic one we write[|: 

R'kM = ^ \/^^^P(-^^/^) - ■ (6) 

The indices "a" and "6" refer to the minimum and maximum of the potential V{Q), located 
at Qa and Qb, respectively. The Ef, stands for the height of the barrier Ej, = V{Qb) — V{Qa)- 
The factor ^jMa/Mb, not contained in Kramers' original work, is meant to account for the 
modification one gets for variable inertia. Notice, that this inertia both influences the current 
over the barrier as well as the number of "particles" (phase space points) sitting in the well. 
Commonly both quantities are calculated with the same M which then drops out; see e.g. 

^ Since we are only looking at stationary situations we leave out the time dependent factor which sometimes 
is taken into account to simulate the "transient time" it takes before the stationary current has built up. 
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Figure 1: Kramers' correction factor to the rate as function of r/: the sohd and dotted curves 
correspond to the high and low (H) viscosity hmit (the latter is calculated for a the barrier 
height of El) = 5.8 MeV, for ^^^Th at 1 MeV) and the dashed ones to the overdamped limit 



eqs. (4.30) and (4.31) of |T0|]. General reasons for the presence of this additional factor will 
be given in Appendix^. At first, in sect . [A. 1| we follow the more common derivation involving 
the construction of the densities at the barrier and at the minimum, which in sect. |A.3| is 
reduced to the Smoluchowsky limit (see below). In sect. |A.4| we follow the arguments of 
Strutinsky [14] which lead to exactly the same formula. Unfortunately, in [H] this feature 
is disguised by the very fact that in the final expression eq.(16), like in large parts of the 
derivation, primed and unprimed quantities are interchanged. It is for this reason that we feel 
compelled to redo the short calculation. The form (^) has recently been derived also in [|1^ 
by applying a generalized version of the "Perturbed Static Path Approximation (PSPA)". 

Notice please that the dependence of the rate on the effective damping strength rjf, still 
is given by (|^). Using the relation ro^ = \C\/M (see (|^)) the limiting value at zero damping 
Rx^'iVb = 0) may be written in the two equivalent forms 



■K 



2^ 



^exp(-i../T) 




exp{-E,/T). 



(7) 



The influence of dissipation is visualized by plotting in Fig.0 the ratio R^^' (Vb) / Rr"' iVb = 0). 
In addition to the result of formula we also show two other cases. Firstly, we show a 
simplified version of the low viscosity limit, 



R^K^ivb 



Rrivb = o) 



It is valid for very small viscosity only (see eq.(|TT|) below) and provided |]T0| the action for the 
motion on top of the barrier may be approximated by /(-Eb) ~ Eb/wb- As demonstrated in 
[^, for the nuclear case such a situation is found only at very small temperatures, much below 
the critical temperature for pair correlations to become important. Secondly, we explicitly 
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indicate the limit the ratio (|]) takes on for overdamped motion 



f^r^^^^ , = — for r^, » 1 . (9) 

Whenever the "high viscosity hmit" apphes the influence of dissipation manifests itself in 
a reduction of the decay rate over the value given by (0). This deviation is claimed to allow 
for deducing a possible temperature dependence of dissipation through the "measurable" 
rate. It must be noted, however, that for overdamped motion it is not the effective damping 
factor rjh which one deduces. Indeed, overdamped motion is governed by the Smoluchowski 
equation in which no inertia appears (see Appendix ^) . But the latter not only is present in 
rjf, but in R^^\rii, = 0) as well. A better way of writing the rate formula in this case is 



i?ovd = R'^K^'iVb > 1) = Tf Ji^^eM-Et/T) ^ — J^—eM-Et/T) . (10) 






\Cb\ 


\Cb\ 


lb 
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\Cb\ 
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Here, the time scale appears which is relevant for overdamped motion across the barrier, 
see below. As can be inferred from Fig.|I] this limit is actually given for values of r]h just 
above unity. Notice, please, that it is only with the additional factor ^/WjMb included in 
(|^), on top of Kramers' classic version, that the inertia drops out in the overdamped limit. 

A few comments are in order on the validity of the rate formulas of the high viscosity 
limit, for which the following assumptions must hold true: 

• On the way from the minimum to the barrier the temperature must not change. 

• The barrier must be sufficiently pronounced, first of all in the sense that its height be 
large as compared to the temperature, viz Ei, ^ T, for further details see sect.|A.l. 



• The effective damping rate must not be too small 

.*>^. (11) 

otherwise formula (||) would have to be applied. 

It may be quite a delicate matter to fix or calculate the temperature which is at stake here. 
For instance, a temperature Ton associated to the total available energy for the compound 
nucleus might be much larger, as for high initial thermal excitations the system may cool 
down by emissions of neutrons or 7's before it fissions. Finally, we should like to remark 
once more that presently any possible quantum features are discarded, which might show up 
at low temperatures 



3 Microscopic transport coefficients 

Evidently, the temperature dependence of the rate will greatly be infiuenced by that of the 
transport coefficients — on top of the infiuence through the Arrhenius factor exp{—Eiy{T)/T). 
Let us first look into results obtained applying linear response theory within the locally 
harmonic approximation 0, before we turn to discuss other forms used in the literature. 
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like in 0]. In this theoretical approach the transport coefficients of average motion are 
obtained by relating in the low frequency regime the strength distribution of a microscopically 
calculated response function Xggi^) to the one of the damped oscillator. The latter is defined 
as 

{Xosc{uj)y^ q{uj) = - {Mio^ + t^LO - C) q{uj) = -gext(tu) • (12) 

and thus may be obtained by adding to (|l]) the term —qext{t) on the right and performing a 
Fourier transformation. For overdamped motion the response function turns into 

i 1 

Xovd(w) = r-TTT ■ (^3) 

7 07 + 26/7 

In accord with the remarks from above on the Smoluchowski limit no inertia appears any- 
more. 

This approach permits one to calculate the transport coefficients as functions of shape 
and temperature for any given nucleus. The formulation is done in such a way that on top 
of shell effects and pairing (see with references to previous works) coUisional damping 
is accounted for as well (for a review see 0). As one may imagine, such computations are 
quite involved, last not least because much knowledge is required about various aspects of 
the dynamics of complex nuclear systems. This is one of the reasons why as yet numerical 



computations have been done only for particular nuclei or for more schematic cases [16 - 
ig, ||. Nevertheless, this experience may allow us to deduce some gross features which 
may be considered generic to a wider class of nuclear systems. This is what we are going 
to do below. It seems appropriate, however, to ffist add some general remarks concerning 
calculations based on the deformed shell model as an approximation to the general mean 
field. 

The output of calculations of the type just mentioned contains much more detailed in- 
formation than at present one may possibly relate to observable quantities. The coordinate 
dependence of the transport coefficients, for instance, is one prime example. Often in nuclear 
transport theories one simply has aimed at constant coefficients for inertia and friction. If 
calculated within the linear response approach, on the other hand, sizable variations with 
shape are seen. One may recall that a similar feature is already seen in the potential land- 
scape, when calculated with the Strutinsky procedure, for instance. Besides the maxima and 
minima which are typical for gross shell effects one sees detailed fine structure. Such features 
may depend on peculiarities of the underlying shell model, and may thus be unphysical in 
nature already by that reason. For the dynamic transport coefficients themselves further 
implications arise from quasi-crossings of levels. To large extent such effects can be expected 
to become much weaker in a multidimensional treatment, which at present still is infeasible. 

One should not forget that problems of this type are intimately related to the fact that 
transport coefficients of inertia, stiffness and friction are those of average motion, calculated 
on the level of the mean field. Finally, however, they are needed for an equation of motion of 
Fokker-Planck type which accounts for dynamical fluctuations. The latter will help to smooth 
out the variations of the transport coefficients in most natural way. Evidently, the problem 
at stake here reflects the general deficiency of mean field theory. In a more appropriate 
treatment one would be able to treat self-consistently both the mean field as well as its 
fluctuations. Since such a theory is not available we suggest some other, more pragmatic 
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procedure. As described previously already, see e.g. P], one may smooth static energies 
as well as the other transport coefficients with respect to their dependence on deformation. 
The averaging interval in Q is to be chosen large enough to wash out the rapid fluctuations 
but small enough to preserve gross shell structures. 



3.1 Shell effects on potential landscape 

The static energy is calculated in the usual way as sum of a liquid drop part and the 
shell correction, both of which depend on temperature, for details see An example of 
the deformation dependence of this potential energy is shown in Fig.0. The dotted curve 
represents the case of zero thermal excitation. On top of the typical gross shell structure 
fluctuations of smaller scale are recognized as well. Features of this type lead to the rapid 
variations of the transport coefficients we talked about above; for the local stiffness as the 
second derivative of the static energy this is immediately evident. As mentioned, we consider 
such fluctuations as unphysical, for which reason we like to remove them by averaging over 
an appropriate interval AQ but to keep the gross shell structure. For the free energy, for 
instance, the smoothing can be done in the following way: 

(14) 

The smoothing function fav{x) in (0) is taken to be that of the Strutinsky shell correction 
method and the Qi are some points in deformation space. The use of the Strutinsky smooth- 
ing function guarantees stability of the averaging procedure: the smooth component of the 
deformation energy is restored after smoothing again. This implies that the liquid drop part 
of the energy is unchanged by this averaging. 

In the figure we also show the averaged potential corresponding to temperatures T = 
0, 1,2 and 3 MeV. As to be expected, with increasing temperature the deformation energy 
becomes much smoother and the height of the fission barrier gets reduced. This is due to 
the reduction of shell effects, as well as to the temperature dependence of the liquid drop 
part. At temperatures above T ^ 3 MeV the shell effects have disappeared completely and 
the averaged deformation energy coincides with its liquid drop component. As seen from the 
figure, at smaller temperatures the shell correction, albeit averaged, does contribute to the 
deformation energy and, hence, to the stiffness. For example, at T = 1 MeV the stiffness 
at the barrier (maximum of the deformation energy) is still several times larger than that of 
the liquid drop part. 

It should be mentioned that in a somewhat different (averaging) procedure was used. 
There the deformation energy was approximated by two parabolas and the stiffness (at the 
minimum and the barrier) was defined by the curvature of these parabola. In this way shell 
effects are washed out to larger extent, not only with respect to the fine structure but even 
with respect to gross shell features. Consequently, the stiffness defined in this way is rather 
close to the liquid drop stiffness. 



3.2 Temperature dependence of transport coefficients 
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Figure 2: The deformation dependence of the collective potential energy. The dotted curve 
shows the deformation energy at T = 0. The dashed, dot-dashed, solid and solid with 
stars curves correspond to the averaged deformation energy at temperatures T = 0, 1, 2 and 
3 MeV. The deformation parameter here is the distance R12 between the centers of mass 
of left and right parts of nucleus (devided by the diameter 2Ro of the sphere with equal 
volume). The averaging is carried out on the interval A{Ri2/Ro) = 0.1 . 




3.2.1 Local stiffness 

In the following we identify stiffness as the one corresponding to the free energy. Then we 
may write 

C(T) = Cldm + SC{T) (15) 

Here, 5C{T) represents the contribution from the shell correction, which disappears with 
increasing temperatures. To parameterize the latter feature we take over a formula of 
to get 

SCm = SC(T = 0)^— (16) 
smh T 

with the shell correction parameter 

9 r , 41 MeV , , 

r=2.'— aud M2„ = -^ (17) 

being the average shell spacing. Above Tgheii — 3 — 5, which corresponds to a temperature of 
the order of 



- shell 



(3_5):^ ~ l-2MeV^ 



277 



the SC{T) practically vanishes, such that C{T) attains its liquid drop value. Eventually, 
this Cldm may still be treated as T-dependent |Q. Finally, we should recall our suggestion 
from above to average the transport coefficients over smaller intervals in Q. In this sense 
the 5C{T = 0) is meant to only represent gross shell features. 
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3.2.2 Local inertia 



As mentioned already in the Introduction, one should expect the inertia to vary with tem- 
perature. It is more than tempting to assume a form similar to the one for the stiffness, 
namely 

M(T) = Mldm + SM{T) (19) 



in which the last term drops to zero like given in (piq). Indeed, within the linear response 
approach a behavior of that type has been observed in a numerical study [|I^]. There, the 
value reached at larger temperatures was given by that of irrotational flow, which for the 
present notation means to identify Mldm = Mrrot- To the best of our knowledge, there is no 
other theoretical model where such a transition is seen explicitly — although one must say 
that in phenomenological applications of transport models commonly the M^j-j-ot is taken to 
represent the macroscopic value of inertia. At present the conjecture behind (p!9D still lacks a 
direct and general proof. However, in the nucleonic response function has been studied 
applying Periodic Orbit Theory (POT). There it was seen that its "fluctuating part" 5x(^) 
decreases with T like the shell correction to the static energy. For slow collective motion 
the inertia is determined by the second derivative of this response function with respect to 
frequency calculated at u; = (see App.^ and eg. ( p.8|) in particular). Therefore, within such 
a model the "shell correction" to the inertia was indeed proven to behave as claimed above, 
although several questions remain open. Amongst others, it is unclear to which extent 
this proof would get modified after considering "collisional damping". The latter cannot 
be treated within POT but should play a major role for the transition to hydrodynamic 
behavior. Possible reasons for rendering a microscopic approach quite difficult have been 



reported in [|1^, [T9|, g. On the microscopic level they are related to the strength distribution 
for the local collective motion. The liquid drop model, on the other hand, represents motion 
of a system having a sharp surface, in contrast to microscopic calculations involving the 



diffuse surface of the mean field and, hence, of the density, for details see |19|. Fortunately, 



however, at larger T when the collisions become more and more important the motion gets 



strongly damped such that the inertia drops out anyway, see (p^Sf). 

Because of these difficulties in microscopic computations we propose to fix the M(T) 
through the vibrational frequency w and the local stiffness by the relation given in (^, 
namely M = \C\/w'^. For our version of Kramers' rate formula this was easily achieved in 
using the second variant shown in (0). 

3.2.3 Vibrational frequency 

At the extremal points of the potential landscape this frequency g7 is a well defined quantity. 
To some extent it is even accessible to experimental verification, at least for zero thermal 
excitation. At the minimum it may be associated with the energy of a collective mode (for 
very recent work on this subject see ||2^) and for the barrier it infiuences the penetrability, 
as encountered for instance in neutron induced fission ||2^. Generally, the hw is believed to 



be of the order of 1 MeV. Indeed, numerical calculations for ^^^Th show this to be 

quite insensitive to temperature; to lesser extent this is true also for the variation with shape 
and mass number. Altogether, for a first orientation the following choice seems appropriate 

hwa ~ ntUfe ^ 1 MeV (20) 
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Figure 3: The inverse relaxation time l/rkm = Fkm/^ (lei^ panel) and (right panel) as 
function of temperature: the microscopic results (solid curves) compared to the approxima- 
tions ( pUD and ( prp (dotted curves). 



with deviations being within a factor of 2 or less. This appears to be the case even when 
pairing is included at smaller T. In Fig.R| we take up the case of ^^^Th, again. The calculation 



is the same as reported in 0; more details will be given below in sect. p.2.6| . From the right 
panel it is seen that this conjecture is pretty much fulfilled. 

3.2.4 Ratio of friction to inertia 

As said above, see eq.(§), the ratio 7/M determines the inverse relaxation time to the 
Maxwell distribution. For underdamped motion this quantity also defines the width Fkin 
of the strength distribution. In Fig.^ we show it on the left hand panel as function of T. 
The dashed curve represents the following approximation, details of which are discussed in 
App.0, namely 

-^h = Fkin ^ 2Fs„(u, T) = , , , , ^ -— MeV (T in MeV) . (21) 

As expected it represents the microscopic result quite well at smaller values of T which 
correspond to smaller values of damping. Recall, please, that the overdamped limit is given 
already for values of the damping factor r^f, slightly above 1, see Fig.|^.. As can be inferred 
from Fig.^ and the estimate (^), this happens at temperatures above 2 MeV; mind that 
77=(7/M)(2a;)-i. 

3.2.5 Ratio of friction to stiffness 

In Fig.^we plot, as function of T, the time Tcou = 7/|C*|, which measures the local relaxation 
in the coordinate. We may recall from (^ that for the overdamped case this is the only 
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Figure 4: The relaxation time Tcou for collective motion at the potential minimum and at the 
barrier: the microscopic result (solid curve) compared to the approximation (|2^ - p^ (dotted 
curve). 



relevant time scale left. Its inverse determines the width of the strength distribution along 
the imaginary axis (see (p!3D). Likewise the decay rate ([T0| ) associated to the Smoluchowski 
equation is proportional to r~i\. In Fig.|, again, the fully drawn line shows the microscopic 
result. The dashed curves represent an approximation, into which the following two features 
are incorporated, the decrease of the stiffness (as given in (|15D and (|16D) and the fact that 
with increasing T the friction coefficient reaches a plateau |18[. To combine both effects 
we chose a functional form similar to the one for the Fkin of (^TD (see App.|BD but with a 
different cut-off parameter Cmacro: 

-JC\" ^ 1 + ^^TV4_ " 1 + vr^TVcLcro^ ^ ' ' 

(22) 

with ~ 1 MeV. One should expect the Tcoh to reach a macroscopic limit like 



|C(r)| 



7waii/2 ^23) 



\ChDM(T)\ 



at larger temperatures. With a parameterization as in (|22D the limit is obtained above 
Th.T — Cmacro/TT, for which reasou the Cmacm would be given by 

2 ^ hw^Tp 7„aii/2 _ 7waii MeV^ 

2 |Cldm(T)| ^ ■ |Cldm(T)| n ■ ^ ^ 



Here, we accounted for results obtained by several previous numerical calculations, see e.g. 



p!9| , g . They showed that the value of friction at large T is somewhat below the wall formula. 
The factor 1/2 is only to be considered a rough rule of thumb. For the stiffness, on the other 
hand, the macroscopic limit evidently is given by the liquid drop model. As the microscopic 
calculation was done with a T-dependent |Cldm(7')| we chose the same one in this fit. In 
both curves the effects of pairing were included, which we are going to address now. 
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Figure 5: The influence of pairing on tlie inverse relaxation time Ty_-^^ = Fkin/^. The micro- 
scopic results of are compared with approximation (pSl) . 



3.2.6 The influence of pairing 

This problem has recently been studied in [Q. The fully drawn lines shown in the previous 
figures refer to such a calculation. Whereas in |0 one concentrated on the regime in which 
pairing is expected to be effective, the present results extend up to T = 4 MeV. Calculations 
in that regime had been reported before in |T^. The underlying shell model is the same in 



both cases, but a different procedure is applied for the single particle width Fgp. For the 
unpaired case the form ( |B.11|) was used, for which the frequency dependence of Fsp(co',T) 



leads to convolution integrals in the response functions. They are known to reduce the 
collective widths 0. In the paired case such a calculation is no longer feasible, for which 
reason there a constant Fsp(/i, A,T) had been assumed, with A being the pairing gap. To 
have a more or less smooth transition to the unpaired case, we now approximate the Fsp(ci;, T) 
by the Fsp(/i, A, T) which above the critical temperature for pairing reduces to the Fsp(/i, T) 
given in ( [B.13| ). For this reason our present friction coefficient may be overestimated slightly. 
For Fkin/^ = 7/M the new results are shown in Fig.^ where we concentrate on temperatures 
up to 1 MeV. To simulate the apparent effect of pairing to reduce friction we suggest the 
modified formulas 

j_ _ 7(A = 0) ^ , rHn(A = 0) .... 
M ~ ^P"" M(A = 0) " ^P"" h ^ ' 

and 

1^ = Tkin = /pair T-km(A = 0) . (26) 

Here, /pair parameterizes the decrease of friction due to pair correlations: An ansatz like 

/pair — ; 7 7^ rr \\ ' (27) 

l + exp(-a(T-To)) 

may do with the following parameters: a = 10 MeV^^ and Tq = 0.55 MeV at the barrier, 
and a = 12 MeV~^, Tq = 0.48 MeV at the minimum. It was found that this choice fits best 
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the microscopic results (with the functional form (^5|) ). It is seen that these values vary with 
shape, which may be an indication that they are perhaps different for different nuclei. For 
a first orientation such details might be discarded. Then an average of the two values could 
be used both for a and Tq. Finally, we should like to remark that this influence of pairing is 
most dramatic for friction, and much less so for M and C . Therefore, we suggest to neglect 
the influence on the latter. 



4 Temperature dependent decay rates 

If one wants to gain information on the transport coefficients and their T-dependence, in 
particular, one needs to separate the influence of the pre-factors from the more or less trivial 
exponential part (the "Arrhenius factor"). Commonly, this feature is then simply identified 
by Kramers' conventional factor which only depends on A systematic study performed 
in |]2B[ has revealed the appearance of a threshold temperature Tthresh above which deviations 
from the statistical model are seen over a wide range of fissioning systems. Moreover, its 
ratio over the temperature dependent barrier heights i?Bar(^) = Eh{T) showed a remarkable 
insensitivity on mass number A. For later purpose it is more convenient to divide that ratio 
by two to get 

T 



~ 0.13 from Fig. 4 of pg. (2^ 

thresh 



We are now going to view this result in the light of our microscopic transport coefficients. 
At first we shall concentrate on the ratio B!^^'{r]i,) / R^^\rji, = 0), to comment later on the 
T-dependence of the additional factors seen in (^, and which involve ratios of the inertias 
or stiffnesses. 



4.1 Rate from microscopic transport coefficients 

In Fig.p we plot the normalized rates for the three cases shown already in Fig.|I], but now as 
function of temperature as determined by our microscopic ri{T). Various conclusions may 
be drawn from this figure: 

• Above T ~ 0.6 — 0.7 MeV an "onset of dissipation" is seen, indeed, in the sense of 
a decrease of the rate for the "high viscosity" limit. The value of this "threshold 
temperature" is in reasonable agreement with the Tthresh — 1-09 — 1.22 given in Tab.l 
of 11 for 224Th. 



Moreover, from the right panel of Fig.^ together with the value for the frequency given 
in Fig.^ one may deduce for these temperatures the r^b to be larger than the value of 
T/2Eb given in (p8[) . It then follows from ( pT]) that the deviation from the statistical 



model is to be attributed to the high viscosity limit. 

This feature evidently is related to the disappearance of pair correlations at T ^ 
0.5 MeV. The curve for the "low viscosity limit" demonstrates that below this tem- 
perature the former case should not be used at all (mind condition (pT])). 
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T (MeV) 



Figure 6: Temperature dependence through a microscopic rj = 1]{T) (fully drawn line); the 
dashed curve represents the overdamped case as given by (||), the dotted one corresponds to 
the low viscosity limit as given by @ for the same barrier height as in Fig.(|l]). 



As indicated earlier, above T 
suggested to use formula (p!0|). 



2 MeV the overdamped limit applies, for which we 



4.2 Comparison with phenomenological models 

In Fig.^we compare the result shown in the previous figure (for formula (^) with the one of 
a macroscopic picture, together with that suggested in 0, §]. By the macroscopic model we 
mean to use the wall formula for friction, the liquid drop model for inertia (irrotational flow) 
as well as for the stiffness. In the following functional form for //^(T) had been suggested: 



r]{T) = 0.2 + 3T^ 



(29) 



It is seen that for such a T-dependence the onset of dissipation occurs quite abruptly and at 
rather small temperatures. One must have in mind, though, that the ratio plotted in Fig.|^ 
is sensitive to the temperature at the barrier, as it is this one which determines the rib- This 
fact might explain why the transition exhibited here through the dashed curve lies at smaller 
T than discussed in p, |^. It can be said that for the other ansatz ri(T) = 0.2 + 5T used 
as well the transition would occur at even smaller temperatures. The reason seems obvious: 
In these two forms the dissipation rate rj is much too large at small T — at least too large 
as compared to our microscopic results. Indeed, as discussed above, at small temperatures 
pairing correlations require dissipation to vanish. 



4.3 Relation to the statistical model 

The rate of the high viscosity limit when extrapolated down to zero friction, R^^irn, = 0), 
corresponds to a variant of the transition state method, and thus to one of the Bohr- 
Wheeler formula (see also fl^)- We demonstrate in sect. |A.2| that this still holds true 
for the modified version of variable inertia. One must have in mind, however, that such an 
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T (MeV) 



Figure 7: Temperature dependence of the decay rate in high viscosity hmit: Fully drawn 
curve: same as in Fig.^; dashed curve: rjb given by (|2^) ; dashed-dotted curve: wall friction 
and liquid drop values for stiffness and inertia. 



extrapolation may be meaningless as, by its very construction, the transition state method 
bases on the assumption of a complete equilibrium. The latter is given only if sizable 
friction forces provide sufficiently fast relaxation to this global equilibrium. Our result for 
the temperature dependence of friction then imply that any version of the transition state 
result must be taken with reservation when applied at small thermal excitations. At this 
point it would be too much to compare complicated evaluations of the Bohr- Wheeler formula 
with estimates of Kramers' rate which have our microscopic transport coefficients as input. 
This must be subject of further studies. 

What is feasible, however, is to compare our results with those of a statistical model, in 



which the Bohr- Wheeler formula |27] is approximated by 



Rstat = 7^ exp{-Eb/T) , (30) 



see e.g. [^, This approximate form comes up when in the Bohr- Wheeler formula p^l the 
level density at the barrier is identified as that of the total system, whereas in the correct 
expression the collective degree of freedom has to be excluded (see e.g.[|l3). That such an 
approximation may lead to erroneous results when interpreting data has been pointed out 



also in |^0[. We may briefly follow up this discussion by using our microscopic input. 

The ratio Rk/ Rstat may be calculated for the three cases we looked at before, low and 
high viscosity limit as well as the overdamped limit. For "high viscosity" one gets from (P) 
(mind (0)): 



R^^^ / Ca 



Rstat T \i \Cb 

For overdamped motion, expression ([TD|) leads to 



Rk 

Rstat 



/a 




h 




76 


T 



i + nl-vA. (31) 



(32) 
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Figure 8: The ratio (p2D of Kramers' rate to that of the statistical model (fully drawn line). 
The dotted curve demonstrates the influence of the stiffnesses, which now are calculated 
from the liquid drop model, but with friction unchanged. 



The result ( pOD can be expected to deviate sizably from Kramers' one. This is demonstrated 
in Fig-H by the fully drawn line which represents the ratio given in (0). Here, the "onset of 
dissipation" seemingly is even more pronounced, and the deviation starts at a much higher 
temperature. However, these effects are not only related to dissipation. Besides the more or 
less obvious fact of the prefactor in ( pO| ) being not identical to the frequency HzUa — hzuh, there 
appears the square root of the ratios of the two stiffnesses. The latter is largely influenced by 
shell effects, which are known to be sensitive to variations of temperature. The implication 
of this feature on the ratio Rx/Rstat is exhibited in Fig.|| by the dotted curve. Its deviation 
from the solid one is solely due to the stiffnesses being evaluated from the (T- independent) 
liquid drop model Cldm(7' = 0). 

Obviously, the Rstat becomes small at small temperatures, say below about 0.5 MeV. In 
this range nuclear friction is small, too, not only at the barrier but also inside the well (see 
Figs.^ and ^). Discarding any quantum effects, which in this regime may become important 
[^, one might then use Kramers' low viscosity limit. As mentioned previously, any transition 
state result, however, does not apply, for which reason a comparison of both is meaningless. 
For the Bohr- Wheeler formula to be valid the system inside the well has to be in complete 



equilibrium [10|. Such a situation is given only for sufficiently large damping. 



4.4 Isotopic effects 

In different nuclei have been studied experimentally and analyzed with respect to a 
temperature dependence of the dissipation strength rjh (called 7 there). In particular two 
isotopes of Thorium have been examined, namely ^^^Th and ^^^Th corresponding to neutron 
numbers of = 134 and 126, respectively. The different behavior seen in experiment has 
been fully attributed to this r]f, and in this way very different values of rjf, have been found, 
together with a different increase with T. Evidently such features cannot be explained within 
a macroscopic picture. One needs to account properly for shell effects. From our experience 
with microscopic computations it seems unlikely that they would have such a big influence 
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Figure 9: Results for different Th isotopes; left panel: Ratio of Kramers' rate to that of the 
statistical model, right panel: influence of shell structure through ratio of stiffnesses. 



on friction alone. On the other hand, recalling that = 126 corresponds to a closed shell, 
it is evident that the stiffnesses for N = 134 and 126 will be quite different, in particular for 
the ground state minimum. This effect can be seen on the right panel of Fig.^, which shows 
the square root of the relevant ratios to be quite different for the two isotopes. Moreover, it 
is easy to convince oneself that this effect is in the right direction. Indeed, as can be seen 
from (pT]), smaller values of this ratio ^/C^JJCb\ simulate larger values of rji,. Notice, please, 
that an excitation of about 100 MeV, as given in Fig.9 of 0, corresponds to a temperature 
of T 2 MeV. Finally, on the left panel in Fig.^ we show the full ratio as given by 
Evidently, we are not able to fully explain the results shown in Fig.9 of 0, but this would 
have been asking too much, for various more or less obvious reasons. 



5 Summary and discussion 

One of our main goals was to suggest simple forms of the temperature dependence of those 
quantities which parameterize transport in collective phase space. These suggestions are 
based on simplified pictures for the intrinsic response combined with experience in micro- 
scopic computations within realistic approaches employing the deformed shell model. Rather 
than using the transport coefficients themselves we argued in favor of combinations which 
allow for a more direct physical interpretation. 

These results were tested at the fission decay rate. To make our reasoning as transparent 
as possible, a simple one dimensional model was used, for which Kramers' picture can be 
applied. This model is somewhat schematic, in so far as the potential is assumed to have 
one pronounced minimum and one pronounced barrier only, which follow more or less closely 
the form of two oscillators, one upright and one inverted. Different to the original work of 
Kramers and of many subsequent applications, our transport coefficients were allowed to 
vary along the fission path. In a first approximation this leads to a modification of Kramers' 
original rate formula. For variable inertia, such a modification already appears necessary for 
the simple reason that for large damping one should reach the rate formula corresponding to 
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the Smoluchowski equation. Still, the deviation from the undamped case is solely determined 
by the damping strength r]h calculated at the barrier. However, for a complete understanding 
of temperature effects of the decay rate it does not suffice to concentrate only on this damping 
strength. 

Of course, such a simple model for the static energy may not apply at all at (smaller) 
temperatures when shell effects lead to important deviations. It is unclear how one could 
generalize the formulas of Kramers and Bohr- Wheeler to realistic cases with two or perhaps 
three barriers having minima in between. To calculate the decay rate for such potentials it is 
perhaps simpler to use the more general description in terms of Fokker-Planck or Langevin 
equations. It may then also be possible to account better for the variations of the other 
transport coefficients with coordinate and temperature, namely inertia and friction. To the 
best of our knowledge so far in such calculations only the transport coefficients of macroscopic 
models have been used. However, the content of Fig.^ indicates that greater deviations can 
be expected for microscopic inputs. 

Nevertheless, within our model we are able to clarify a few important points which must 
not be discarded even at temperatures at which shell effects do not really dominate the 
process. 

• The transition to overdamped motion already occurs at temperatures around T ~ 
2 MeV. Then the decay rate should be calculated from the Smoluchowski equation. 

• In this case no inertia appears any more. Thus neither the frequency nor the effective 
width Fkin = h'-f/M (sometimes referred to as hjS) play any role. 

• The solely important quantities are the ratio of friction to stiffness, 7/C, and (for the 
decay rate) the ratio Ca/Cb of the stiffnesses. Physically, the former determines the 
relaxation time for sliding motion in the potential, see (|22D. In passing we may note 
that it is essentially this time which determines the saddle to scission time, an effect 
not considered here. Notice that these stiffnesses are subject to large shell effects which 
may easily be accounted for by the shell correction method. As a first demonstration 
of this effect we have calculated the impact of these stiffnesses on the decay rate of the 
Thorium isotopes N = 134 and = 126. 

• In principle, the vibrational frequency w should be obtained from microscopic com- 
putations (for T < 2 MeV). But a simple and quite fair estimate may be given by 
hzua — hzub — 1 MeV (|20|) , which is almost exact for the case of ^^''Th. Evidently, the 
zu may be greatly influenced by shell effects, but from our experience we claim that 
the deviations are within a factor of about two. 

• Finally, we should like to briefly comment on the paper [^. There scaling rules have 
been derived for the case of the Smoluchowski equation, and based on phenomenological 
input. In the schematic model underlying the discussion a few assumptions had been 
made which are not in accord with microscopic results. First of all, the ratio of barrier 
height to temperature definitely decreases with T. Secondly, as the authors monitor 
themselves, a frequency at the barrier of the order of 20 MeV is much too high. Together 
with the value used for P, namely /3 = T]^i^/h ~ 10^^ sec~^ one gets an rjb of ~ 0.2. This 
implies the motion around the barrier to be under damped rather than overdamped. 
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Moreover, for several cases studied in this paper, the rj}, is seen to be of the order of 



or even smaller than T/2Eh. Thus, according to ( jTT]) not even Kramers' high viscosity 
limit seems appropriate. 

Let us turn now to very low temperatures, say to the regime where pairing correlations 
become important. They imply an additional reduction of dissipation, such that one may 
truly speak of the onset of dissipation when pair correlations disappear. It must be said 
though that this transition occurs at temperatures definitely smaller than those suggested in 
[m " Hi- mentioned earlier, one has to make sure, however, that one speaks of the same 
temperature. The one of the compound nucleus might be larger than the one the system 
still has when it passes the barrier. 

It ought to be stressed that a small damping strength at small temperatures may have 
quite drastic implications. In case the "high viscosity limit" still applies quantum corrections 
to the decay rate would lead to an increase of the latter Q]. If, however, the dissipation 
strength falls below the limit given by ([Tl|) the nature of the diffusion process would change 
completely. Then dissipation is too weak to warrant relaxation to a quasi-equilibrium. This 
not only violates Kramers' rate formula (for the high viscosity limit), or our extension of it, 
but also the Bohr- Wheeler formula becomes inapplicable. Moreover, so far no method exists 
of how one might incorporate collective quantum effects. 

Finally we should like to indicate that our transport coefficients are not free of uncer- 
tainties in some of the parameters specifying our microscopic input. The most difficult one 
is found in the single particle width. It has been parameterized in by referring to the 
optical model for nucleons in finite nuclei through a kind of local density approximation. 
There the Tq and its dependence on the nuclear density had been traced back to microscopic 
computations of the self energies in a generalized Brueckner type description. This implies 
an uncertainty of perhaps a factor of two in all formulas where the friction coefficient ap- 
pears, as in (^) or (22) as well as in the r] of (|). Since any microscopic answer on such 



questions is extremely difficult, there is hope of narrowing down this uncertainty by more 
elaborate comparisons with experimental results. 



A Implications from a variable inertia 



In Kramers' seminal paper [0 the equation for the density in phase space was written and 
applied to the decay of a metastable system for the case of constant transport coefficients. 
In nuclear physics both the inertia as well as the friction coefficient vary with the collective 
coordinate. This has been accounted for already in early applications to heavy ion collisions, 
where globally Gaussian solutions (centered at the classical trajectories) were used to cal- 
culate reaction cross sections see e.g. [Q. Written in compact form for one variable the 
transport equation should read: 



d_ 
di 



f{Q,P,t) 



+ 



MiQ) dQ dQ \2M{Q) + dP 



(A.l) 



f{Q,P,t). 
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Neglecting any quantum effects, which might show up at small temperatures only |^ , the 
diffusion coefficient is given by the classic Einstein relation D{Q) = 'y{Q)T. A possible T- 
dependence of the transport coefficients has not been indicated explicitly. This is immaterial 
as long as we treat temperature as a fixed parameter, which is assumed to hold true in the 
entire paper. As is easily verified, one stationary solution of ( [A.l|) is the distribution of 
global equilibrium, 

/eq(Q, P) = ^ exp{-pn{Q, P)) (A.2) 
with the energy given by the classical Hamilton function 

The reason is due to the following features: (i) The conservative part of the equation has 
the form of Liouville's equation, namely 

- f{Q, p, t) = {niQ, p) , /(g, p, t) } . (A.4) 

(ii) The terms which represent dissipative and fluctuating forces are assumed independent 
of the momentum P, such that P only appears quadratically in the kinetic energy, (iii) The 
choice of the diffusion coefficient by D{Q) = 'y{Q)T makes the second line of ( |A.1| ) vanish 
once applied to (|A.2|) . 

The structure of these equations not only allow for the proper equilibrium, it also warrants 
the continuity equation to be valid in the form 

g-^n{Q,t) + —j{Q,t) = 0, (A.5) 

with the current and spatial densities being defined as 

j{Q,t)= I dPjJ^^f{Q,P,t) (A.6) 

and 

n{Q,t)= [ dPf{Q,P,t), (A.7) 



respectively. This is easily verified with the help of ( |A.1| ) exploiting partial integrations with 



respect to momentum P (for which no "surface terms" survive). 

Both in Liouville's equation as well as in the transport equation ([A.l| ) a term appears 



which needs to be treated with special care. It is the one which results from the spatial 
derivative of the kinetic energy, and reads 



Indeed, this term is absent in the common derivation of the rate formula (see e.g. [|10|), 
where always a constant inertia is assumed to be given. It may be noted that this term 
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may imply difficulties with "saddle point approximation", as needed in Kramers' stationary 
solution. Moreover, and more important for the present purpose, this term is also neglected 
in extracting the transport coefficients within the linear response approach. There, a locally 
harmonic approximation is exploited, which for the sake of simplicity is formulated only with 
respect to the coordinate. What that means may be visualized in the following way. Look 
at the Liouville part of the transport equation, in particular at the term which involves the 
derivative of the Hamilton function with respect to Q. First of all there is the ordinary force 
from the potential, which in the expansion around a Qq to harmonic order may be written 
as 

--^ViQ) ^ K{Q) ^ K{Qo) - C{Qo){Q - Qo) , (A.9) 

with the second derivative of the potential defining the local stiffness C{Qo). In addition 
there is the term 

dQ 2M{Q)) dQ \2M{Q)) ^ ' ^ 

which is of second order (namely in P). In a consistent treatment one would have to introduce 
a Pq and expand all terms locally in collective phase space around the (Qo, Pq) to second 
order both yhQ — Qq and P — Pq (see ^3[)- this has not been done when extrapolating 

transport coefficients from the microscopic linear response approach. This does not imply 
that our inertia may not change with the coordinate at all. It only means that the terms 
including its derivative have been discarded. With respect to the basic transport equation 



( |A.1| ) this approximation may be phrased as 

d 



d P^ \ d 
dQ2M{Q)) aP^^^'^'^) 



(A.ll) 



A.l Kramers' decay rate formula extended to variable inertia 



For reasons given above, we still want to make use of the condition ( |A.11| ). Nevertheless, it 



is necessary to re-derive an expression for the rate, which still is defined as 

i? = A. (A.12) 

To calculate the quantities involved we need a global solution fg\oh{Q,P) of equation (|A.1|) 
which corresponds to a small but finite, (quasi-) stationary current across the barrier. This 
current may be calculated from (|A.6|) , with f{Q,P,t) replaced by fgiohiQb, P)- The proba- 



bility Na of finding the system inside the well at Qa may be calculated as follows 

Qa+A Qa+A 

Na= J dQ J dPfg,oh{Q,P)= J rfQngiob(g). (A.13) 

Qa-A Qa~A 

The integration range 2A has to be smaller than Qh — Qa but large enough such that it 
contains the vast majority of the ensemble points sitting in the well. An approximation for 



22 



fg\oh{Q,P) may be constructed by matching together at some intermediate point Q local 
solutions valid at the minimum Qa and at the barrier Qj,. The global solution fgioh{Q,P) 
will have an overall normalization factor which drops out when calculating the rate from 
(A.12). The normalization of the local solutions f{Q ~ Qa,P) and f{Q ~ Qb,P), on the 



other hand, might be different. It has to be chosen in appropriate fashion such that both 
solutions match properly at Q. 

For a sufficiently high barrier the particles inside the potential well may be assumed to 
stay close to a local equilibrium associated to the temperature T. In case that the corre- 
sponding fluctuations {AQ'^)'^ concentrate on a region around the minimum, the potential 
may be replaced by a harmonic oscillator and the associated local phase-space density may 
be approximated by: 

/„(g, P) ^ f{Q ^ Q„, P) = Ma exp (^-/? + V{Qa) + Y (Q - Qaf^ ) • (A.14) 

This is a reasonable estimate up to a Q with 

(Q-Qa)'>(AQ2)r = ^. (A.15) 

Likewise, approximating the barrier by an inverted oscillator, the phase-space density may 
be represented by Kramers' stationary solution 0, (neglecting any quantum effects, see 
e.g. i) 



h{Q, P) = f{Q ^ Qt, P) =A4 exp (^-(3 + v{Q,) - ^ (g - Q,)' 



X I du , exp I I 



(A.16) 



where 



I a 



A = '—^ ^ and a = TM^ ^ - 1 | (A.17) 



This latter solution must be joined to the one given by ( A.14 ) in some intermediate region 



say at g. It would be too much to require this to be possible for any P, in particular as we 
are not able to treat the inertia in continuous fashion, for reasons given above. However, as 
the rate is determined by the ratio of two quantities which are averaged over momentum it 
turns out sufficient to match only the reduced g-space densities. Like suggested by ( [A.7| ), 



the latter are obtained by integrating out the momentum P in ( |A.14 ) and ( A.16|) to get 



riaiQ) = Afa V'^nMaT exp (v{Qa) + - ^-)')) (A.18) 



and 




UbiQ) = Mb \l exp {-p (vm - ^(g - Qbf] ) 1 1 + erf I ^/ ^ (g - g, 

(A.19) 
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To obtain the last expression identities for error functions have been apphed. Still it is not 
of a form for which a condition like na{Q) = nb{Q) would make much sense. For this we 
need the additional assumption 

|a|(Q-Qfe)'>T, (A.20) 

which renders the error-function in ( |A.19| ) close to unity. Then the two densities match 
smoothly at a Q satisfying 



\Cr, 



(Q-QbY 



if only the normalization constants are chosen according to 



(A.21) 



(A.22) 



Now we are in the position to calculate the rate from (|A.12|) . Plugging (|A.22|) into (|A. 14|) 
the number of "particles" at the minimum ( [A.13| ) becomes: 



'a ~ 



Mb^2TiMbT exp( 



Mb 



Afb2nT^I-^ exp( 



-pV{Qa)) (A.23) 



To get this simple expression it was assumed that the A, which in ( |A.13| ) defines the range 
of integration, is of the order of or larger than the fluctuation {AQ'^)'^ given in ( |A.15D , such 
that the Gaussian integral can be calculated for A ^ oo. The current at the barrier may 
be evaluated from ( [A.16|) and ( [A.6|) with Q = Qb- After a lengthy but straight forward 
calculation involving identities for error-integrals, once more, one arrives at 



Jb 



AfbT 



Mb_ 

\Cb\ 



VI + Til - Vb exp{-pV{Qb)) 



(A.24) 



From (pr^ and the decay rate { \KJ^ turns into 



Tjh.v. " " 




Tjb exp{-pEb) 



Tjb I exp(-/5E{, 



(A.25) 



confirming the expression (^) used in the text. We may note again that the result ([OsD 
was derived before in within an extension of the Perturbed Static Path Approximation 
(PSPA). 

Finally, we like to come back once more to the conditions ( [A. ISj ) and ( |A.2U| ) imposed 
before. They go along with the relation ( [A.21| ) for Q and the barrier height Eb = V{Qb) — 
V{Qa). The latter must then satisfy Eb = ^(Q - QaY + ^(Q - Qbf > T. It may be 
useful to visualize these relations with the help of the following schematic potential: 

V(Qa) + f (g-Qa)' for Q<g, 

(A. 26) 



V{Q) 



Eb + V{Qa) 



\Cb 



{Q-Qb? for Q>Q. 



Choosing the Q according to Q = {CaQa + \Cb\Qb) / {Ca + \Cb\) the two parabolas match 
smoothly with a continuous first derivative. Possible errors related to (|A.15|) and (|A.20|) 
may easily be estimated from elementary properties of the error function. 
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A. 2 The relation to the transition state result 



In transition state theory one assumes to be given a system that is totally equilibrized 
inside the barrier and for which at the barrier current only flows outward discarding any 
backflow. Within its most general version, the fission rate has been estimated by N. Bohr 
and J. A. Wheeler through their famous formula ||2^. There the equilibrium is the one of a 
micro canonical ensemble as represented by the density of states (of the total system at the 
minimum and of the intrinsic system at the barrier). To the extent that the micro canonical 
ensemble may be represented by a canonical one, with one and the same temperature at the 
minimum and at the barrier, the calculation of the rate can be done as follows, looking only 
at the collective degree of freedom. The outward current (at the barrier) is given by 



■trans 
Jb 



poo p poo p / p2 \ 



Comparing with Kramers's stationary solution shown in ( |A.16| ) one realizes the only 
difference being the replacement of the Theta function 0(P) of (|A.27|) by the integral (for 
Q = Qb) which appears in the second line of ( |A.16| ). Because of the following representation 
of the Theta function 



Q{P) = hm / du -= exp -— , (A.28) 
"^^J-oo V27rcr V 2,a J 

it is seen that (for finite temperature) 

Rtrans=RK''-{'n = ^)- (A.29) 

This follows immediately with the help of the expression given for a in (A. 17). 



A. 3 The Smoluchowski limit 

Performing in ( |A.25|) the limit r/b ^ 1, with the effective damping rate defined by (^, one 
gets: 

S /Si "^"^-^^'^ = i (A.30) 

Rovd • 

This expression coincides with the formula ([lOD associated above with the Smoluchowski 
limit. As a matter of fact, this result can be obtained directly from the Smoluchowski 
equation 

d_ g / 1 dV{Q) T d 

The calculation of the decay rate is quite easy in this case. Indeed, eq. ( |A.31| ) gives an explicit 
form of the j{Q,t) as a function of the density n{Q,t). Although in our case the friction 
coefficient varies with Q the common derivation of the rate formula (see e.g. [^) may be 
taken over without much difficulties. 



(A.31) 
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Notice, please, that the transition from an equation hke (|A. 1|) to (|A.31|) only requires the 
r]{Q) to be large enough at any Q. Such a transition may be performed also in the case of a 
variable inertia, at least if condition ( |A.11[ ) is fulfilled. In any case, in the overdamped limit 
the inertia has to drop out. We may note in passing that this transition is in accord with 
the locally harmonic approximation in the form discussed in sect. 2. 2. 5 of 0. Following the 
arguments of sect. 10.1 and 10.4 of |^ eq. ( |A.31| ) can strictly be derived from ( [A.l| ) neglecting 
the term {\K^ ). 



A. 4 Strutinsky's derivation of the rate formula 



The essential idea exploited in |jT^ is written there below eq.(6). Different to the approach 
described in sect. |A.l|, the number of particles Na at the minimum is estimated by multiplying 
the density nb{Q) of Kramers' stationary solution calculated at Qa by an "effective length", 
which in turn is determined by the mean fluctuation of the oscillator at this Qa times v^27r, 
viz by 



cq 



2'kT 



(A.32) 



(The additional factor \^2tt is required to ensure the appropriate measure needed for the 
normalization of a Gaussian). In this way one gets from and ( |Al9| ) 



jyh.v. 




vl~Vb] exp{-(3Eb) 



(A.33) 



Here, it was assumed (i) that the Qt — Qa is sufficiently large such that in ( |A.19| ) the 
error function could be replaced by unity, and (ii) that the barrier height can be estimated 
as Eb ^ {\Cb\/2){Qa - QbY- The result (|A.33|) has the same form as given in (|A.25|) . As 
explained earlier, it is equivalent to (^ or the second line of (|A.25| ) which involve the inertias. 



This latter expression is identical to the one given in eq.(16) of |]T4| if one only interchanges 
there primed and unprimed quantities. 



B A schematic microscopic model 
B.l The Lorentz-model for intrinsic motion 

Let us assume that the nucleonic excitations can be parameterized by the response function 
(in this section we set h = 1) 



^ [u;-n + iT/2 u + n + iT/2 

Here, the average matrix element of the one body operator F, which acts as the generator 
of collective motion, measures the overall strength of the distribution. The states reached 
by that coupling are centered at fl with an effective bandwidth F (measured here in units 
of MeV/h). For real frequencies the reactive and dissipative response functions, x' and 



(B.l) 



26 



x", are readily calculated noticing that they represent real and imaginary parts of x{^) = 
x'{uj) + ix"{uj). The static response is given by 

It is useful to rewrite this (intrinsic) response function x{^) terms of the form of the 
oscillator response given in (|12]). 

^^'^^ + zFo; - + (r/2)2) uj^ + tr,^tOO-r:ol,- 

In this way transport coefficients for intrinsic motion appear 

Mi,t = -^, ri„t = r, zul, = n' + {T/2f. (b.4) 

Zilr ^ 

Next we turn to the collective response. For the F-mode it is given by ||^ 

X{uj) 1 -l/Mj. 

1 + kx[uj) ^ + k uj^ + iTpuj - Wp 

with the inverse coupling constant 

-i = C(0) + x(^ = 0); (B.6) 

and C(0) being the stiffness of the free energy. The transport coefficients for the collective 
i^-mode are: 

M^ = Mi,t, r^ = ri,t, Cf = MfwI = M,^,wl, + k (B.7) 

To get those for the Q-mode one needs to multiply these quantities by For slow modes 

it so turns out that to a good approximation the C(0) in ( [B.6|) may be neglected as compared 
to the x(0)- This leads to 



(B.8) 

a;=0 



k^ 2QF^ (fi2 + (r/2)2)' n^ + {T/2y 2 doo^ 

In the last expression we have made use of (|B.2|). Last not least this has been done because 



the static response seems to be quite insensitive to the increase of temperature |^6[, at least 
for not too large T. The transformation from the F- to the Q- mode leaves ratios between 
transport coefficients unchanged. The collective width, the ratio between friction and inertia, 
thus becomes Fkin = F. For friction this implies 



F dx" 
7 = FkinM = ^^^rn^^ x(0) 



+ (F/2)2 ' duj 



(B.9) 

oj=0 
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and for the stiffness one gets the expected result C ~ C(0), which follows because the Cp of 
( |B.7| ) can be written as 



X(0) x(0)(c(0)+x(0)) (;^(0)) 



The second equation follows from (B^) and 

Finally, we should like to note that for the schematic model with only one mode the 
inertia always is the one which defines the value of the energy weighted sum. Likewise, as 
one may see from (p.9| ) for friction and (to lesser extent) from ( p.8| ) for the inertia, these 
transport coefficients are well represented by their "zero frequency limits". 



B.2 Benefits and shortcomings of this model 
B.2.1 Weak damping 

For this model weak damping is defined as ^ r/2. In this case static response and inertia 



turn into the expressions known from the so called "degenerate model" x(0) = 2F'^ /VL 
and M = 2F'^/Q^. Notice that the inertia has the typical structure of the cranking inertia. 
For friction one gets 7 = 2rF'^/Q^. 

The degenerate model becomes most transparent if it is applied to the case that nucleons 
move in oscillator potentials, in particular if any spin dependent forces are neglected. Then 
the intrinsic excitation is given by hn = ANhQo = AN (41 MeV/A^^^), where AA^ is the 
difference in the major quantum numbers of those states which are coupled through the 
multipole operator F. Whereas for the quadrupole there is only one possibility, namely 
AA^ = 2, this is no longer true for other multipoles, for which more than just one mode are 
possible. The same holds true as soon as a spin orbit force is introduced. Then even for the 
quadrupole transitions with AA^ = are possible. It is them which lead to the low frequency 
modes we are typically interested in, as they resemble closest the fission mode. If one still 
likes to stick to the (degenerate or) Lorentz model — which only allows for one mode — the 
effective frequency Q will only be a fraction of the shell spacing parameter Qq. 



B.2. 2 Strong damping 

It is tempting to apply this schematic model also to the extreme case of very strong damp- 
ing when F becomes comparable to or larger than the frequency Q of the typical intrinsic 
excitation. Plain confidence into the formula (p.9| ) would lead to 7 ~ (4/F)x(0). This 
seems particularly intriguing if on the one hand the static response does indeed not change 
much with T, and if, on the other hand, the F is associated to the widths the single particle 
states, as will be discussed below. According to ( |B.ltj| ) there might then be some range in 
which the friction force would show the typical dependence one expects for liquids in 
the "collision dominated regime", see also sect. 5. 3 of Jl^. However, we claim that for finite 
nuclei the situation is more complicated. Evidently, the effects of strong collisions are due to 
the increasing importance of residual interactions. But the latter imply other consequences 
as well, last and not least a mixing with more complicated states such that with increasing 
thermal excitations many particle - many hole states become more and more important. As 
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has been demonstrated in previous papers, see e.g. [|T3], || amongst others, this effect imphes 
that high frequency modes shift to lower frequencies such that the typical mode at stake in 
the transport model gets more and more strength — implying that finally its inertia is given 
by the sum rule limit. Moreover, it has been demonstrated that this feature goes along with 
the disappearance of shell effects at T = Tsheii. This problem is addressed in the text. 



B.2.3 Temperature dependence through collisional damping 

Looking back at the intrinsic response function introduced in ( |B.1D , one realizes that the 
only quantity which can be expected to change sensitively with excitation is the width F. To 
get some first orientation we may relate it to the single particle width. For a Fermi system 
the latter can be expected to be of the form H] F| 



with the parameters 



= 0.03 MeV"^ and c = 20 MeV . (B.12) 



For slow collective motion we may omit the frequency dependence and evaluate this width 
at the Fermi surface and put — /i = in (p.ll| ) . Along this approximation we may put 



0.6T2 

1 + TV40 



Tkin ~ 2rsp(/i, T) ^ ^ 'rr.9.;.r. ^eV (T in MeV) . (B.13) 



Evidently the correction term in the denominator only becomes important at temperatures 
of the order of T ~ 6 MeV. This is already beyond that value were the other effects come 
into play we discussed in sect . [B. 2.^ . For this reason the actual r(T) is changed in the main 
text. 

Finally we may note that our schematic model is not capable of accounting for pairing. 
The latter will modify the transport properties at temperatures below T ~ Tpair- This is 
discussed in the text. 
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